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Abstract. We study the Log(N)-Log(S) and X-ray luminosity function in the 2-10 keV energy band, and the 
spatial (3-D) distribution of bright, Lx > 10 34 — 10 35 erg s _1 , X-ray binaries in the Milky Way. In agreement 
with theoretical expectations and earlier results we found significant differences between the spatial distributions 
of low (LMXB) and high (HMXB) mass X-ray binaries. The volume density of LMXB sources peaks strongly 
at the Galactic Bulge whereas HMXBs tend to avoid the inner ~ 3 — 4 kpc of the Galaxy. In addition HMXBs 
are more concentrated towards the Galactic Plane (scale heights of w 150 and « 410 pc for HMXB and LMXB 
correspondingly) and show clear signatures of the spiral structure in their spatial distribution. The Log(N)- 
Log(S) distributions and the X-ray luminosity functions are also noticeably different. LMXB sources have a 
flatter Log(N)-Log(S) distribution and luminosity function. The integrated 2-10 keV luminosities of all X-ray 
binaries in the Galaxy, averaged over 1996-2000, are ~ 2 - 3 • 10 39 (LMXB) and ~ 2 - 3 • 10 38 (HMXB) erg 
s" 1 . Normalised to the stellar mass and the star formation rate, respectively, these correspond to ~ 5 • 10 28 erg 
s" 1 Mq 1 for LMXBs and ~ 5 • 10 37 erg s _1 /(M yr" 1 ) for HMXBs. Due to the shallow slopes of the luminosity 
functions the integrated emission of X-ray binaries is dominated by the ~ 5-10 most luminous sources which 
determine the appearance of the Milky Way in the standard X-ray band for an outside observer. In particular 
variability of individual sources or an outburst of a bright transient source can increase the integrated luminosity 
of the Milky Way by as much as a factor of ~ 2. Although the average LMXB luminosity function shows a break 
near the Eddington luminosity for a 1.4 Mq neutron star, at least 12 sources showed episodes of super-Eddington 
luminosity during ASM observations. We provide the maps of distribution of X-ray binaries in the Milky Way in 
various projections, which can be compared to images of nearby galaxies taken by CHANDRA and XMM-Newton. 
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1. Introduction 

Recently the CHANDRA X-ray observatory studied the 
distributions and luminosity functions of X-ray binaries 
in at least 7 spiral, e. g M 81 (Tennant et al. 2001 ), 2 ellip- 
tical, e.g. NGC 4697 flSarazin et al. 2000| ), and 2 starburst 
galaxies, M 82 ( Zezas et al. 2001 ) and Antennae ( Fabbiano 
et al. 200 1[ ). The main discovery of these CHANDRA ob- 



servations was the existence of numerous point-like sources 
with luminosities in the CHANDRA spectral band consid- 
erably higher than the Eddington luminosity of a 1.4 M 
neutron star. Nearby galaxies observed by CHANDRA 



Send offprint requests to: Hans- Jakob Grimm, e-mail: 
gr imm@mpa-gar ching . mpg . de 

Correspondence to: Max-Planck-Institut fur Astrophysik, Karl- 
Schwarzschild-Str. 1, 85748 Garching, Germany 



have a great advantage compared to observations of X-ray 
sources in our Galaxy: All objects observed in a particular 
galaxy are equidistant and therefore it is straightforward 
to construct the luminosity function in the CHANDRA 
band. However, even with the angular resolution and sen- 
sitivity of CHANDRA we are restricted to nearby galaxies 
(d & 50 Mpc) and we are able to observe only the high 
luminosity end of the luminosity function. 



Observations of compact sources inside our Galaxy 
thus open the unique possibility to construct a luminosity 
function in a much broader range of luminosities and this 
might be important to construct the synthesised spectrum 
of the LMXB and HMXB populations of the Galaxy in a 
broad spectral range from 0.1-500 keV using data from all 
existing spacecraft. 
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In this paper we use data of the All-Sky Monitor 
(ASM) ( [Levinc et al. 1996|) ab oard the Rossi X-ray Timing 
Explorer (Brandt et al. 1996) to investigate the following 
topics. 

— Using ASM data, existing information about the 
source distances and a model of the mass distribution 
in the Milky Way we constructed the luminosity func- 
tion of high and low mass X-ray binaries in our Galaxy. 

— Distribution and number of high mass X-ray binaries 
are expected to trace the location and reflect the rate 
of star formation. The X-ray luminosity of starburst 
galaxies might become an additional source of infor- 
mation about the star formation rate in these galax- 
ies. Moreover knowledge of the luminosity of HMXBs 
versus star formation rate opens the way to find the 
volume cmissivity of our universe at different redshifts 
in hard X-rays due to starburst and young galaxies. 

— The luminosity of the LMXB component is propor- 
tional to the total mass of the old stellar population of 
the Milky Way. Concerning the volume emissivity of 
galaxies these data will provide information about the 
contribution of elliptical galaxies and old star popula- 
tions in spiral galaxies only at sufficiently low redshifts 
(z < 0.4 — 0.5). To obtain the volume cmissivity due 



than the Eddington value for a neutron star. In at 
least 7 sources the compact object has been iden- 
tified as a neutron star based on the detection of 
X-ray pulsations or X-ray bursts, therefore we know 
with certainty that the peak luminosity exceeded the 
Eddington limit. Moreover, the total number of super- 
Eddington sources might be higher because we know 
from the broad band observations that the bulk of the 
luminosity can be emitted outside the ASM sensitivity 
band. 

In terms of the spatial distribution of X-ray binaries 
this paper elaborates on works done earlier that also dis- 
tinguished between low and high mass systems but used 
substantially smaller samples. 



Previously White et al. (1980), Lamb et al. (1980), 



Nagase (1989) and Vcrbunt (1996) noted the correlation 
of the positions of accreting X-ray pulsars with high mass 
companions with the location of spiral arm features of the 
Milky Way. Based on a larger sample of HMXBs with 
measured distances we show that indeed the spatial dis- 
tribution of HMXBs follows the spiral structure of the 
Galaxy. 

Using distance estimates and angular distribution of 



var 



to q>ld star populations at higher redshift we need to 



LMXBs |van Paradijs fc White (1995[ ) and [White 
Paradijs (1996| ) investigated the spatial distribution of 



know a model of the mass exchange rate evolution. 
Our Galaxy should become an important point in 
the future calibration curves of Lhmxb/SFR and 

-^LMXB/-Wgalaxy 

Our analysis of ASM data permits us to show how 
our Galaxy would look from outside in different pro- 
jections. This will allow us to compare data about our 
Galaxy with new CHANDRA observations. 
Surprisingly enough, just the comparatively few most 
luminous Galactic X-ray binaries practically dominate 
the X-ray luminosity of our Galaxy. The majority of 
the brightest X-ray binaries are extremely variable on 
all time scales from milliseconds to years-tens of years. 
Therefore the luminosity of our Galaxy as a whole 
would also be subject to strong variability. This is 
important because with a powerful X-ray telescope 
such as XEUS it will be possible to detect X-ray flux 
from distant galaxies on the level of L ~ 10 40 erg 
s _1 but only short time scale variability would permit 
to distinguish the collective emission of X-ray binaries 
from the low luminosity, AGN-type activity of the nu- 
cleus. Black holes are unable to produce strong vari- 
ability with characteristic times significantly shorter 
than a few 0.01s ^ flSunyaev fc Revnivtsev 2000] ). 
For super-massive black holes the characteristic time 
is of order or above ~ 10 3 s. 

Our analysis of ASM data and data from other space- 
craft shows that at least for 17 X-ray sources in our 
Galaxy ASM or other spacecraft detected flux reaching 
or exceeding the level corresponding to the Eddington 
critical luminosity for a 1.4 Mq neutron star, see Table 
|l|. Maximal fluxes detected were up to 10 times higher 



LMXBs and BHC in our Galaxy, particularly in the 
Galactic disk. They estimated values for the vertical (290 
pc and 710 pc for BHC and NS binaries) and radial scales 
(4.5 kpc for NS binaries) of the disk. These values are 
in general agreement with those obtained in this paper, 
that are based on a considerably larger number of sources. 
jGrcbenev et al. (1996 ) found good agreement between 
the source distribution observed by ART-P/GRANAT in 
the Galactic Centre region and the stellar mass distribu- 
tion in the Galactic Bulge. We thus have a reasonably 
good knowledge about the distribution of LMXBs in the 
Galaxy. 

2. RXTE All-Sky Monitor Data 

In order to construct the Log(N)-Log(S) distributions and 
luminosity functions we used the publicly available data 
of ASM. The ASM instrument is sensitive in the 2-10 keV 
energy band which is divided into 3 broad energy channels 
and provides 80% sky coverage for every satellite orbit (~ 
90 minutes). Due to its all-sky nature and long operational 
time, ~ 5 years, the ASM instrument is ideally suited for 
studying time averaged properties of sources. The light 
curves are obtained by RXTE GOF (Lcvine et al. 1996) 
for a preselected set of sources from the ASM catalogue. 
The catalogue consists of sources which have reached an 
intensity of more than 5 mCrab at any time (Lochner & 
Rcmillard 1997), and as of June 2000 included 340 sources 
of which 217 are galactic and 112 extragalactic, and 10 
unidentified. The distribution of galactic sources on the 
sky is shown in Fig. [l|. For a detailed description of selec- 
tion criteria and a list of sources see Lochner & Rcmillard 
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Table 1. Persistent, transient and extragalactic X-ray binaries for which episodes of near- or super-Eddington flux were 
detected. The maximum luminosities from ASM observations refer to the dwell-time light curves (90 s observations 
every 90 minutes) . Thus the values might differ from the luminosities given in Tables ^| and |^. 



Source 


type 


Mi 


Luminosity [10 38 erg s 1 ] 


Energy range 


RefW 


distance 


Ref. (c) 






[M ] 


average 


a) peak 


[keV] 




[kpc] 




Persistent sources 


Cir X-l 


LMXB 


NS 


4.4 


12 


2-10 


(1) 


10.9 


(0 


GRS 1915+105 


LMXB 


14-30 


3.7 


15 


2-10 


(1) 


12.5 


(ii) 


Sco X-l 


LMXB 


NS 


2.7 


9.4 


2-10 


(1) 


2.8 


(iii) 


Cyg X-2 


LMXB 


NS 


1.8 


4.2 


2-10 


(1) 


11.3 


(i),(iv),(v),(vi),(vii) 


GX 349+2 


LMXB 


NS 


1.6 


3.2 


2-10 


(1) 


9.2 


(i),( v ii),( vi ii) 


GX 17+2 


LMXB 


NS 


1.5 


3.0 


2-10 


(1) 


9.5 


(i),(vii),(ix),(x) 


GX 5-1 


LMXB 


NS 


1.4 


2.2 


2-10 


(1) 


7.2 


(i),(vii) 


GX 340+0 


LMXB 


NS 


1.3 


2.2 


2-10 


(1) 


11.0 


(i),(vii) 


Cyg X-3 


HMXB 


NS(?) 


0.5 


2.1 


2-10 


(1) 


9.0 


(xi) 


X 1624-490 


LMXB 


NS 


0.24 


3.3 


2-10 


(1) 


13.5 


(ix) 


GRO J 1744-28 


LMXB 


NS 


0.15 


4 


8-20 


(2) 


8.5 


(xii) 


Transient sources 


V4641 Sgr 


HMXB 


9.6 




33 


2-10 


(1) 


9.9 


(xiii) 


GS 2023+338 


LMXB 


12 




11 


1-40 


(3) 


4.3 


(xiv) 


4U 1608-52 


LMXB 


NS 




9.2 


2-20 


(4) 


4.0 


(i),(x),(xv) 


N Muse 91 


LMXB 


7 




6.1 


1-6 


(5) 


5.5 


(xvi) , (xvii) , (xviii) 


XTE J1550-564 


LMXB 


10.5 




5.3 


2-10 


(1) 


5.3 


(xix) 


N Oph 77 


LMXB 


5 




5.4 


2-18 


(6) 


7.0 


(xviii), (xx) 


GS 2000+251 


LMXB 


6 




2.2 


1-6 


(7) 


2.7 


(xviii), (xxi) 


Magellanic Clouds sources 


SMC X-l 


HMXB 


NS 


2.0 


17 


2-10 


(1) 


60^ 




LMC X-l 


HMXB 


4.7\/cos i 


1.5 


13 


2-10 


(1) 


50^ 




LMC X-2 


LMXB 


NS 


1.5 


17 


2-10 


(1) 


50< d > 




LMC X-3 


HMXB 


> 5.8 


1.5 


17 


2-10 


(1) 


50< d > 




LMC X-4 


HMXB 


NS 


0.38 


15 


2-10 


(1) 


50< d > 





average luminosity observed by ASM. 
^ Reference for the peak luminosity. 
(1) ASM (this paper ), (2) 

(1978|), (7) frsunemi et al. 



Watson et al 
1 Rcferencc(s) for the distanc e 



Sazonov et al. (1997| ), (3) |Tanaka (199^ ), (4) |Nakamura et al. (1989| ), (5) [Kitamoto et al. (1992j ), (6) 



(i) Ivar 
(viPff 



1989 



Mirabel . 



Rodrig 

(iv) |Orosz fc Kuulkers (1999|), (v) pawley et al. (1979[), (vinSmale (1998|), (vii) |Pc nniii\ 



(i) |van Paradi.js fc White (1995|), (ii 



(1996y (ix) phristian fc Swank (1997^ 
Orosz et al. (2000[), (xiv) |King (1993f). 



(x) [Ebisuzaki et al. (1984|) , (xi) |Predehl et al. (2000[) (xii 
tv) |Nakamura et al. (198£ ), (xv i) Preincr et al. (1994) 



(19941) 



(iii) 



Bradshaw et al 



(1999|) 



(198S), (viii) |Wachter fc Margon 



Barret et al. (1996[ ), (xix) prosz et al. (2002| ), (xx) |Martin et al. (1995[ ), (xxi) phevalier 
w assuming a distance of 50 kpc for LMC and 60 kpc for SMC. 



Nishiuchi et al . (199S| ), (xiii) 
(19960 , ( xviii ) 



ii) Orosz et al 
Ilovaisky (1990| j 



(1997). The 1 day sensitivity of ASM is w 10 mCrab corre- 
sponding to a count rate of 0.75 cnts s _1 . The ASM count 
rate has been converted to energy flux assuming a Crab- 
like spectrum and using the observed Crab count rate: 

F[erg s _1 cm" 2 ] = 3.2 • 10~ 10 • i?[cnts s -1 ]. (1) 

The 1-dwell ASM light curves have been retrieved from the 
RXTE public archive Q at HEASARC and cover a time pe- 
riod from the start of the mission through 27/04/00. In or- 
der to construct Log(N)-Log(S) the light curves have been 
averaged over the entire period of available data which 
might differ for different sources. We did not account in 
any way for orbital variations or eclipses, as e.g. in Cen 
X-3. 

1 ftp:/ /legacy.gsfc. nasa.gov/xte/data/archive/ASMProducts/ 
definitive.! dwell/ 



Important for the analysis presented below are the 
questions of systematic errors in the light curves and of 
the completeness limit of the ASM catalogue. 

2.1. Systematic errors 

The ASM light curves are assumed to have a systematic 
error at the level of ~ 3% which is added in quadrature 
to the statistical errors in the light curves provided by the 
RXTE GOF. The systematic error has been estimated us- 
ing Crab data and refers to the ~ dwell-day time scales. 
The formal errors for the average fluxes calculated from 
the entire ASM light curves are very small ~ 0.1 — 0.2 
mCrab (~ 1 — 2 • 10~ 2 cnts s _1 ). In the presence of sys- 
tematic errors this might not correctly characterise the 
accuracy of the average flux estimate, especially for weak 
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Fig. 1. Distribution of LMXBs (open circles) and HMXBs (filled circles) in the Galaxy. In total 86 LMXBs and 52 
HMXBs are shown. Note the significant concentration of HMXBs towards the Galactic Plane and the clustering of 
LMXBs in the Galactic Bulge. 



Table 2. List of the sources used to estimate systematic 
errors. 



Source 


average fhix (a) excess RMS (6) 
[cnts s _1 ] [cnts s _1 ] 


Cas A 


4.9 ± 0.007 


~ 0.08 


Tycho SNR 


1.3 ± 0.007 


~ 0.04 


Puppis A 


0.84 ± 0.008 


~ 0.05 


Vela pulsar 


0.75 ± 0.008 


~ 0.01 


CTB 33 


0.35 ± 0.014 


~ 0.07 


PSR 1259-63 


0.18 ± 0.012 


~ 0.01 


NGC 2024 


0.09 ± 0.008 


~ 0.02 


PSR J1713+0747 


0.07 ± 0.015 


~ 0.01 


PSR 1957+20 


0.06 ± 0.012 


~ 0.02 


XTE J1906+090 


0.04 ± 0.011 


~ 0.03 



("'the errors are formally calculated using the errors in the 
light curves. 

^ upper limit on the unaccounted contribution of the system- 
atic errors to the averaged flux, estimated from Fig. ^. 



sources. The contribution of systematic errors to the aver- 
age flux estimate depends on their statistical properties, 
in particular their correlation time scale. In order to inves- 
tigate these properties we selected several sources believed 
to have constant X-ray flux, like SNRs or rotation pow- 
ered pulsars, see Table |], and rebinned their light curves 
with different bin durations ranging from 1 to 200 days. 
For each binned light curve we computed the expected 
RMS from the errors given with the light curves and com- 
pared it with the observed RMS. The results are shown 



in Fig. g. Ideally there should be a one-to-one correspon- 
dence between expected and observed RMS (straight line 
in Fig. ||). As can be seen from Fig. || this is not the case. 
The observed RMS somewhat exceeds the expected value, 
the discrepancy increasing towards large bin durations (~ 
50-200 days). The excess variance at large bin durations 
(lower-left part in Fig. ^) gives an upper limit on the un- 
accounted systematic error in the averaged flux estimate. 
As can be seen from Fig. ^ the particular value of the 
systematic error, though varying from source to source, is 
in the range of 0.01-0.1 cnts s _1 . We assumed a value of 
0.05 cnts s _1 (to be added in quadrature to the statisti- 
cal error). We further verified that our conclusions are not 
sensitive to the value of the systematic error. 



For 15 sources we obtain statistically significant, > 
3<7, negative average count rates. The majority of these 
sources, namely 14, are located in the Small and Large 
Magellanic Cloud and their negative average flux is ap- 
parently caused by source interference in these crowded 
regions. The remaining source also appears to suffer from 
interference with nearby sources. In particular, we have 
noticed that some of the light curves show a clear drop 
below zero count rate coincident in time with addition 
of new sources located nearby to the ASM catalogue. All 
these sources are excluded from our analysis. 
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Fig. 5. Number-flux relation for galactic X-ray binaries. The vertical dashed line corresponds to our completeness 
limit of 0.2 cnts s _1 . The solid lines are the best fit models to the ASM data - a power law for HMXBs and a power 
law with cutoff in the differential Log(N)-Log(S) distributions at 110 cnts s _1 for LMXBs (see Eqs.(||) and (||)). 



Table 3. The best fit values for the number-flux relation for different classes of galactic sources from the ASM 
catalogue. 



Subsample 


no. of sources (1) /all sources 


cutoff [cnts s ] 


normalisation 


slope 


quality of fit 
(K-S test) 


all galactic 


131/217 (2) 


110 


88 


0.34 ± 0.05 


92% 




132/217 




72 


0.41 ± 0.04 


51% 


LMXB 


83/105 (2) 


110 


83 


0.2 ±0.06 


71% 




84/105 




56 


0.3 ±0.05 


0.5% 


HMXB 
SNR 


25/51 
6/7 




9.4 
4.8 


o.ei±S:iS 

0.36tg; 22 


46% 
98% 


CV 


5/10 




0.5 


1.68 ±0.61 


98% 



Number of sources above the completeness limit of 0.2 cnts s . 
For fits with a cutoff the brightest source, Sco X-l, was excluded. 



2.2. Completeness 

Important for the analysis presented below are two aspects 
of completeness: 

1. completeness flux limit of the ASM sample of the X- 
ray sources 

2. completeness of the sample of galactic X-ray binaries 
which are optically identified and for which distance 
measurements are available 

The first problem arises for example in studying Log(N)- 
Log(S) distribution of all galactic sources and is addressed 
below. The second problem is important in analysing 
Log(N)-Log(S) distributions of various types of galactic 
X-ray sources and especially their luminosity functions. It 
is discussed in Sect. ||. 

Due to the present method of construction of the ASM 
catalogue its completeness limit is difficult to assess in any 



straightforward way. By definition the ASM sample in- 
cludes all sources, galactic and extragalactic, which have 
reached an intensity of 5 mCrab at any time, which corre- 
sponds to a completeness limit of ~ 0.37 cnts s _1 . On the 
other hand we know from the same ASM light curves that 
non-transient Galactic X-ray binaries have typical values 
of the ratio of maximum flux (on the time scale of dwell-~ 
day) to average flux of the order of few. Therefore, in terms 
of long term average values the ASM catalogue might be 
complete down to lower fluxes. 

In order to indirectly probe the completeness limit 
of the ASM sample we use the fact that the Log(N)- 
Log(S) relation for extragalactic sources is well known and 
follows a power law with index —3/2 (Forman ct al. 1978), 



down to ~ 3.8 • 10" 14 erg s" 1 cm" 2 (Ogasakaet al. 1998) 



which corresponds to ASM count rate of 1.2 • 10 cnts 
s _1 . The Log(N)-Log(S) relation for extragalactic sources 
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A XGC 202-1 
B Vela SNR 

PSR 1259 63 
PSR 1957-20 




A PSR J1713 
□ XTE J 1906 
q CTB 33 

Puppis A 
«m, Tycho SNR 
Cas A 



0747 
090 



expected RMS [cnts/s] 

Fig. 2. Observed versus expected RMS for 10 different 
sources and for different time binnings. The bin duration 
varies from dwell time scale, i.e. ~ 90 seconds (upper right 
corner), to 200 days (lower left corner). Although there is 
considerable spread, the observed RMS is generally higher 
than expected, especially at large bin durations exceed- 
ing 50 days (expected RMS < 0.1 cnts s -1 ). Assuming 



that systematic and statistical errors are independent the 
systematic error may be added to the statistical error in 
quadrature. This is shown by the solid curves for three 
different values of the systematic error: 0.01, 0.05 and 0.1 
cnts s _1 . 



based on ASM data is compared with HE AO A-l and 
ASCA results in Fig. One can see that flattening of 
the source counts caused by incompleteness of the sample 
begins at a count rate of ~ 0.1 cnts s^ 1 . 

Therefore we set, somewhat arbitrarily, the complet- 
ness limit of the ASM sample of the X-ray sources at 0.2 
cnts s _1 . We verified that our conclusions are not sensitive 
to the exact value. 



3. The Log(N)-Log(S) distributions 

In order to calculate the number-flux relations the ASM 
light curves were averaged over the entire time span of 
available data for each source. The resulting Log(N)- 
Log(S) relation for galactic sources is shown in Fig. [|. The 
differentiation between galactic and extragalactic sources 
was done using SIMBAD database. The overall shape and 
normalisation of the Log(N)-Log(S) relation of Galactic 



sources is similar to that obtained by UHURU (Forman 



et al. 1978[ ) and ARIEL V (|Warwick ct al. 1981] ). The 
UHURU result (Matilsky et al. 1973) is schematically 
shown in Fig. [| by the solid line. The Log(N)-Log(S) re- 




flux [cts/s] 

Fig. 3. Log(N)-Log(S) distribution of extragalactic 
sources. Magellanic Cloud sources have been omitted. 
The upper histogram contains all extragalactic sources, 
the lower histogram excludes 4 nearby galaxy clusters 
(Perseus, Virgo/M 87, Coma and Centaurus). The shaded 
region shows the Log(N)-Log(S) obtained by HEAO-1 A- 



2 for high latitude (|6| > 20°) sources (Piccinotti et al 



1982 ). The width of the shaded region roughly accounts 



for the uncertainty of the RXTE/ASM and HEAO-1 A-2 
calibration. 



lation for different types of Galactic sources is also shown 
in Fig. £|. 

We further selected X-ray binaries from the sample 
and divided them into low mass (LMXB) and high mass 
(HMXB) binaries according to the mass of the optical 
companion, using the mass of the secondary, M2, of 2.5 
Mq to separate high and low mass systems. The precise 
value of this boundary affects classification of only few 
X-ray binaries (Her X-l, GX 1+4, GRO J1655-40 etc.). 
In doing so we us ed SIMBAD datab ase, the Catalogue 
of X-ray Binaries ( van Paradijs 1994), the Catalog ue of 
CV, LMXB and related objects ( [Ritter fc Kolb 1998j), the 



catalogues of low-mass X-ray b inaries (Liu et al. 2001 ) 
and high-mass X-ray binaries ( Liu et al. 2000 ) and in 
some cases publications on individual sources. Recently 
the donor star in GRS 1915+105 was identified to be a 



K or M giant (Greiner et al. 2001) so this source is clas- 
sified as an LMXB. Of 115 galactic X-ray binaries with 
average ASM flux exceeding our completeness limit of 
0.2 cnts s _1 only 6 sources were left unclassified. The 
fraction of unclassified sources is ~ 5% and they have 
fluxes in the 3 • 10 x — 13 cnts s^ 1 range and therefore 
should not affect our conclusions in any significant way. 
The compilation of galactic X-ray binaries with type, op- 
tical companion, average flux and, if available, distance 
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Fig. 6. The distribution of Galactic HMXBs (solid lines) and LMXBs (thick grey lines) against Galactic latitude b 11 
(left panel) and longitude I 11 (right panel). The distribution against b 11 of HMXBs shows a stronger concentration 
towards the Galactic plane compared to LMXBs. Along I 11 LMXBs show a strong concentration in the direction 
towards the Galactic centre. The arrows in the right panel mark the positions of the tangential points of spiral arms. 
The broad hump in the HMXB distribution at I 11 — 100° — 160° is mostly composed of relatively low luminosity 
sources in the Perseus and Cygnus arms. Note that on the right panel the number of LMXBs is divided by 3. 



and average luminosity is available in electronic form at 
http://www.mpa-garching.mpg.de/^grimm/. The result- 
ing Log(N)-Log(S) relations for LMXBs and HMXBs are 



shown in Fig 



gradual steepening of the Log(N)-Log(S) relation occurs 
towards higher fluxes. Similar behaviour was also found 
by UHURU QMatilsky et al. 1973| ) and OSO-7 flMarkert 
et al. 1979). We therefore modified the simple power law 



To fit the observed Log(N)-Log(S) distributions we 
used the usual power law in the form: 



JV(> S) = k-S- a 



(2) 



where N(> S) is the number of sources with fluxes higher 
than S, a is the slope, and k the normalisation. S is mea- 
sured in ASM cnts s _1 . In order to calculate the best fit 
values of the parameters we use a Maximum-Likelihood 



method in the form suggested by Murdoch & Crawford 
(1973j)] This implementation of the M-L method takes into 
account the errors associated with the flux. Since the sys- 
tematic error dominates the averaged flux error we used 
the value of 0.05 cnts s _1 from Sect. 

of the error. The error is assumed to be Gaussian. 

-l 



as an estimate 
Only 

sources with an averaged flux above 0.2 cnts s _1 were 
used in the fit. The best fit values for different types of 
Galactic sources are given in Table ^. The errors given 
are an estimate of the la errors for one parameter of in- 
terest derived from the Maximum-Likelihood method. In 
order to characterise the quality of the fit we used the 
Kolmogorov-Smirnov test. 

As is obvious from Fig. |j| and the results of the K-S 
test (Table |^) a simple power law distribution does not de- 
scribe the observed Log(N)-Log(S) relation for LMXBs. A 



in the form: 



N(> S) = k- (S- a - S~ a ax ) 



(3) 



This corresponds to a cutoff in the differential Log(N)- 
Log(S) relation at flux S — S max . The value of the cutoff 
was chosen to S max = 110 cnts s _1 . The results, however, 
are not very sensitive to the actual value of S max . The 
above value of S max corresponds to the ASM flux from 
a 1.4 M Q neutron star located at a distance of 6.5 kpc 
(average distance of LMXBs from the Sun) and radiating 
at Eddington luminosity. For fitting the Log(N)-Log(S) of 
all galactic and LMXB sources with cutoff we excluded 
the brightest source, Sco X-l, from the sample since its 
flux is far higher than the cutoff. As can be seen from 
Table |^ and Fig. || introduction of the cutoff significantly 
improves the quality of the fit for LMXBs. On the other 
hand it does not change significantly the results for other 
types of Galactic sources, especially HMXBs. Note that 
the steepening of the Log(N)-Log(S) for LMXBs is not 
an artifact of the incompleteness of the source sample at 
low fluxes. The numbers do not change qualitatively if we 
increase the low flux limit by a factor of 2 - the values of 
K-S probability are 6% and 68% for a single power law and 
a power law with cutoff in the form of Eq. (0), respectively. 
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Table 4. The parameters of the standard Galaxy model. 
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Fig. 4. Number-flux relation for all galactic sources de- 
rived from the entire ASM sample. The broken solid line 
shows schematically the number-flux relation for the low 
latitude |6| < 20° sources obtained by UHURU flMatilsky 
et al. 1 973 ). The vertical dashed line shows approximate 
completeness limit of the ASM sample. The thick grey his- 
togram shows the Log(N)-Log(S) for all Galactic sources 
observed by ASM. The four lower histograms show the 
contributions of different classes of sources to the total 
galactic Log(N)-Log(S). 



4. Spatial distribution of X-ray binaries 

Progress in the number of distance determinations and 
identifications of secondary stars in X-ray binaries in the 
last decade opens the opportunity to study the 3-D distri- 
bution of XRBs in more detail than was previously possi- 
ble. Notwithstanding the still relatively small number of 
X-ray sources and the sometimes poor accuracy of dis- 
tance determinations it is now possible to compare the 
observed distribution of XRBs with theoretical expecta- 
tions. This is not only interesting in itself but, because 
of the flux limited nature of the ASM sample, knowledge 
of the spatial distribution is required in order to derive 
the luminosity function. Due to the above mentioned un- 
certainties and the flux limitation of the sample it is still 
not possible to unambiguously determine shape and pa- 
rameters of the XRB distribution. We therefore adopted 
an approach in which we use the standard model of the 
stellar mass distribution in the Galaxy as a starting point 
and adjust, whenever possible, its parameters to fit ob- 
served distributions of low and high mass X-ray binaries. 
As the luminosity function depends somewhat on the as- 
sumed spatial distribution, we verify that variations of the 
parameters, which can not be determined from the data 
do not affect derived luminosity functions significantly. 



4.1. Angular distribution of X-ray binaries 

The all-sky map shown in Fig. [I] demonstrates vividly that 
the angular distributions of high and low mass X-ray bi- 
naries over the sky differ significantly. This fact is further 
illustrated by the angular distributions against Galactic 
latitude and longitude shown in Fig. ||. The figures illus- 
trate the well-known fact that HMXBs are strongly con- 
centrated towards the Galactic plane. In addition drastic 
difference in the longitude distributions of HMXBs and 
LMXBs can be noticed, with the latter significantly con- 
centrated towards the Galactic Centre/Bulge and the for- 
mer distributed in clumps approximately coinciding with 
the location of tangential points of the spiral arms, see e.g. 
Englmaier fc Gerhard (1999| ); | imonson (1976[). 



4.2. Source distances and 3-D distribution of X-ray 
binaries 

In order to study the spatial distribution of X-ray bina- 
ries we collected source distances from the literature. We 
found distances for 140 X-ray binaries from the ASM sam- 
ple. For X-ray binaries with an average flux above the 
ASM completeness limit, used for constructing the lu- 
minosity functions in Sect. 0, distances were determined 
for all but 8 sources. In cases when the published dis- 
tance estimates disagree significantly we used the least 
model dependent estimates or their average. For the com- 
pilation of the source distances see bttp: / / www.mpa] 
Igarching. mpg.de/~grimm/ . The spatial distribution of X- 
ray binaries in various projections is shown in Fig. 



4.3. The Galaxy model 

As a starting point in constructing the spatial distribution 
of X-ray binaries we employ the standard three compo- 
nent model of the stellar mass distribution in the Galaxy 



(Bahcall & Soncira 1980), consisting of bulge, disk and 
spheroid. The parameterisation of bulge and disk is taken 



from Dehncn & Binney (199S) and for the spheroid we 
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Fig. 7. Face-on view of the Galaxy - distribution of low 
mass (open squares) and high mass (filled circles) X-ray 
binaries. The origin of the coordinate is at the Galactic 
Centre. The Sun is located at x=0, y=8.5 (marked by the 
pentagon). The thin solid line shows logarithmic 4- armcd 
(m=4) spiral model with pitch angle of 12° (e.g. Vallcc| 
( 1995QJ ■ The thick solid lines show the spiral model of the 
Galaxy based on optical and radio observation of the gi- 
ant HI I region s ( Georgelin fc Georgelin (1976 ), Taylor fc 
Cordes (1993 )). The fact that the majority of sources is 
located at y > is due to the flux limited nature of the 
ASM sample and incompleteness of the optical identifica- 
tions/distance measurements at the large distances from 
the Sun (see discussion in the text). 



take the model of Bahcall fc Soneira (1980| ): 



pBulge — PO, Bulge 1 ( 



''(I 



-) 7 • exp(- 



(4) 



/ Tm r \ z \\ fr\ 

PDisk = Po,Disk ■ exp( ) (5) 

r r d r z 



PSph, 



P0,Sph, 



exp( 



-M#) 1/4 ) 



(6) 



where po,B«Zge, Po,Disk and po,s P here are the normalisa- 
tions, r is the distance in the plane from the galactic cen- 
tre, z is the distance perpendicular to the galactic plane, 
and R is the distance from the galactic centre in spherical 
coordinates. All distances are in kiloparsec. Meaning and 
values for other parameters are given in Table |J. 

In the standard Galaxy model the mass ratios of 
the components are about 2:1:0.3 for disk:bulge:spheroid. 




5 10 15 

distance from GC [kpc] 
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Fig. 8. Radial distributions of high mass (solid histogram) 
and low mass (thick grey histogram) X-ray binaries. The 
projected distance is defined as y/ ' x 1 + y 2 , where x and y 
are Cartesian coordinates in the Galactic plane, see Fig. 
0. Note that the plotted distributions are not corrected for 
the volume of cylindrical shells (oc r). 



These numbers follow from the model using normalisa- 
tions for the disk, po.Disk — 0.05M Q pc~ 3 , and spheroid 
population, p , sphere = 1/500 ■ po^pisk, observed in the 



vicinity of the Sun ( Zombeck 1990) an d a bulge mass of 
about ~ 1.3- 10 10 M a (pwek et al. 1995|) . All these masses 
refer to baryonic mass in the stars. 

All three components of the standard Galaxy model 
were used to construct the spatial distribution of LMXB. 
The spheroid component with appropriately adjusted nor- 
malisation was used to account for the population of glob- 
ular cluster sources. Based on the observed distribution 
and theoretical expectation that HMXBs trace the star 
forming regions in the Galaxy, only the disk component 
was used for the spatial distribution of HMXBs. 

Several parameters, namely vertical scale height of the 
disk and relative normalisation of the spheroid for the 
LMXBs, can be determined directly from our sample of 
X-ray binaries. For these parameters we used the best fit 
values inferred by the data. For the rest of the parameters 
we accepted standard values for the stellar mass distri- 
bution in the Galaxy. The final set of the parameters is 
summarised in Table |]. 

The disk component of the standard Galaxy model was 
modified in order to account for the Galactic spiral struc- 
ture. The description of the spiral arms is based on the 



model of Georgelin & Georgelin (1976) derived from the 
distribution of HII regions. To include it into our Galaxy 



model we used the FORTRAN code provided by Taylor & 
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Cordes (1993). The spiral arms computed in this way are 
shown in Fig. ^ by thick grey lines. This empirical model 
is close but not identical to a 4 arm logarithmic spiral with 
pitch angle of 12° (e.g. Vallee (1995)) shown in Fig. ^ by 
thin solid lines. 

In the following two subsections we discuss spatial dis- 
tribution of HMXBs and LMXBs in more detail. 

4.4. High mass X-ray binaries 

The angular distribution of HMXBs in Fig. || shows sig- 
natures of the Galactic spiral structure. These signatures 
are clearly seen in the distribution of sources over galac- 
tic longitude which shows maxima approximately consis- 
tent with directions towards tangential points of the spiral 
arms. No significant peak in the direction to the Galactic 
centre is present. The signatures of the spiral structure be- 
come more evident in the 3-D distribution of the smaller 
sample of sources for which distance measurements are 
available, Figs. 0J|. The radial distribution (Fig. ||) shows 
pronounced peaks at the locations of the major spiral arms 
and is similar to that of primary tracers of the Galactic 
spiral structure - giant HII regions (e.g. Downes et al 



(1980)) and warm molecular clouds (e.g. Solomon et al 
(1985)). In particular, the central ^3 — 4 kpc region of 



the Galaxy is almost void of HMXB well in accordance 
with the radial distribution of the giant HII regions and 
warm CO clouds. This appears to correspond to the inte- 
rior of the 4-kpc molecular ring. 

The vertical distribution of HMXBs is significantly 
more concentrated towards the Galactic Plane and suffi- 
ciently well described by a simple exponential with a scale 
height of 150 pc as shown in the left panel of Fig. ||. 

Based on theoretical expectations and on the data 
shown in Fig. 0J|^,^ we included only the disk compo- 
nent in the volume density distribution HMXBs. It is clear 
however that a simple exponential disk is not a good de- 
scription for the radial distribution of HMXB. Therefore, 
following Dchnen fc Binncy (1998| ) we assumed the disk 
density distribution in the form given by Eq. (|^), where 
the first term in the exponential allows for the central 
density depression. To describe the observed central de- 
pression for HMXBs a rather large value of r m « 6 — 7 
kpc is required (cf. r m — 4 kpc from Dchnen & Binney 



(1998|)). The spiral arms were assumed to have a Gaussian 



density profile along the Galactic Plane: 



where w a 



PSpirai oc y~] exp(-(^-) 2 ), (7) 
r-i Wa 

600 pc is the width of the spiral arm, and 



Sj is the distance to the nearest point of the spiral arm j 
projected to the Galactic Plane: 



Sj = ^x-x'F + ty-y')*. (8) 

In order to account for the spiral structure the disk den- 
sity, Eq. (||), was multiplied by pspiraV- 

PDisk PDisk ' P Spiral (9) 



4.5. Low mass X-ray binaries 

Contrary to HMXB, the angular distribution of LMXBs is 
strongly peaked in direction to the Galactic centre and de- 
clines gradually along the Galactic plane, see Fig. ^. The 
central ~ 2 kpc region is densely populated with Galactic 
Bulge LMXB sources and contains ~ 1/3 of the LMXBs 
from our flux limited sample (Fig. |J). A noticeable fea- 
ture of the radial distribution of LMXB is the pronounced 
minimum at ~ 3—4 kpc. This minimum approximately co- 
incides with the ~ 1 — 3 kpc gap in the distribution of the 
molecular gas and the ^2.2 kpc minimum in the density 



of infrared light distribution in the Galaxy (Binney et al 



1997) and probably separates bulge sources from the disk 



population. Similar to HMXBs, the signatures of the spi- 
ral structure might be present in the radial distribution 
although they are less pronounced. 

The vertical distribution outside the bulge (Fig. ^) is 
significantly broader than that of HMXBs and includes a 
number of sources at high galactic z. A formal fit to the 
observed distribution with an exponential law results in a 
large scale height of 9 50±130 pc, which is close to th e value 
of 710 pc obtained by |van Paradijs fc White (1995| ) for NS 
LMXBs. However, due to presence of a tail of sources at 
\z\ > 1.5 — 2 kpc, the observed z-distributions cannot be 
adequately described by a simple exponential law. As only 
three out of nine sources &t\z\ > 2 kpc are located in glob- 
ular clusters, this tail of high-z sources cannot be solely 
due to the globular cluster component. A possible mecha- 
nism - a kick received by a compact object during the SN 
explosion, was considered e.g. by van Paradijs & White 
|(1995| ) . The relatively small number of high-z sources does 
not allow one to determine the shape of their distribution 
based on the data only. In order to account for the high- 
z sources and the LMXB sources in globular clusters we 
chose to include in the spatial distribution of LMXBs the 
spheroid component described by a de Vaucouleurs pro- 
file (Eq. (^)). Note that a de Vaucouleurs profile correctly 
represents the distribution of globular clusters. The over- 
all vertical distribution can be adequately represented by 
a sum of an exponential law with a scale height of 410lgQ° 
pc and a de Vaucouleurs profile with the parameters given 
in Table ^[ The spheroid component represented by the de 
Vaucouleurs profile contains a ~ 25% of the total number 
of LMXBs. Note, that this number is by a factor of ~ 2- 
3 larger than the mass fraction of the stellar spheroid in 
the standard Galaxy model. The enhanced fraction of the 
spheroid component is generally consistent with the fact, 
that the number of X-ray sources per unit mass is ~ 100 
times higher in the globular clusters than in the Galactic 
disk and 12 out of 104 LMXBs in our sample are globular 
cluster sources. 

The angular resolution of the ASM instrument does 
not permit to study in detail the very central region of 
the Galaxy which is characterised by the highest vol- 
ume and surface density of X-ray binaries. Based on 
GRANAT/ART-P data having significantly better scn- 
sitivity and angular resolution, Grebenev et al. (1996) 
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Fig. 9. Vertical distributions of high mass (left panel) and low mass (right panel) X-ray binaries. The vertical dis- 
tributions were summed over northern and southern galactic hemispheres. In the case of LMXBs only sources with 
R > 3.5 kpc were used, to exclude bulge sources. The thick grey solid lines show the observed distributions and the 
thin solid and dashed lines the expected distributions for an exponential disk with 150 pc scale height for HMXBs, 
and an exponential with scale height 410 pc and a 25% contribution of the spheroid for LMXBs, respectively. For the 
assumed model see Eqs. (F| ph. 



showed that the distribution of the surface density of X- 
ray binaries in the central 8° x 8° of the Galaxy is con- 
sistent with the stellar mass distribution in the Galactic 
Bulge. 

To conclude, our model of the volume density distribu- 
tion of LMXBs includes all three components of the stan- 
dard model of the Galaxy: bulge, disk and spheroid with 
the disk-to-spheroid mass ratio decreased to 4:1. Similarly 
to HMXBs, r m ss 6 — 7 kpc is required to describe the 
central density suppression of the disk population. The 
modulation of the disk component by the spiral pattern 
at the 20% level was also included: 

PlH^k 3 OC PDisk ■ (1 + 0.2 • pspiral) (10) 

where p Dtsk is given by Eq. (j|) and p Spira i - by Eq. (§). 

4.6. Completeness of the sample of the distance 
measurements. 

The fact that the majority of the sources in Fig. |?] is lo- 
cated at y > is related to the flux limited nature of 
the ASM sample (obviously it is easier to observe weak 
sources located closer to the Sun) and to the incomplete- 
ness of the available distance measurements (more difficult 
to measure the distance to a more distant source). The 3- 
D distribution of X-ray binaries enables one to check the 
latter effect. 



Plotted in Fig. [T(j is the distribution of LMXB sources 
with luminosities Lx > 4 • 10 36 erg s _1 over the distance 
from the Sun. For the ASM completeness flux limit of 0.2 
cnts s _1 , sources with Lx > 4 ■ 10 36 erg s _1 should be 
visible up to a distance of ~ 20 kpc. However, comparison 
with the expected distribution computed using t he L MXB 
volume density distribution constructed in Sect. 4.5 shows 



an increasing deficiency of sources at distances <; 10 — 15 
kpc. In total ~ 14 sources in the distance range of 10- 
20 kpc are "missing". These "missing" sources should be 
hidden among the ~ 20 unclassified sources in the ASM 
catalogue for which no optical identification/distance de- 
terminations are available. 



Recent observations by Kuijken & Rich (2001) lend 



support to this interpretation. They measured proper mo- 
tions of blue and red giants in direction to the Galactic 
centre. The red giants, concentrated in the Galactic bulge, 
have a velocity dispersion in Galactic coordinates, b 11 ver- 
sus I 11 , symmetric around zero. However, blue giants, lo- 
cated in the disk, have a velocity dispersion asymmetric 
around zero with respect to I 11 which means that there 
is a net motion of the observed blue giants in one direc- 
tion. Interpreting this as the motion of the disk around 
the Galactic centre, it also means that there is a deficit of 
the observed blue giants on the far side of the Galaxy (c.f. 
Fig. 0). 

This comparison (Fig. [l0|) shows that our sample of 
optical identifications/distance measurements for LMXB 
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Fig. 10. Distribution of the LMXB sources over distance 
from the Sun (thick grey histogram). Only sources with 
luminosity Lx > 4 • 10 36 erg s _1 are plotted. Given the 
ASM completeness flux limit of 0.2 cnts s -1 , sources with 
Lx > 4 • 10 36 erg s _1 should be visible from the distance 
of upto w 20 kpc. The thin solid histograms shows the ex- 
pected distribution of the sources in the model constructed 
in Sect. 4.5. Deviation of the observed distribution from 
the prediction becomes visible at the distance > 10 — 15 
kpc. 



sources is complete up to a distance of ~ 10 kpc. The 
significantly smaller number of HMXBs above the ASM 
completeness flux limit did not permit us to perform a 
similar analysis for HMXB sources. However, one might 
expect that due to the higher luminosity of the optical 
companion the limiting distance for HMXBs is not smaller 
than for LMXBs. We therefore accepted a value of D max = 
10 kpc as a maximum source distance for the luminosity 
function calculation for both types of sources described in 
the next section. 

5. Luminosity function 

Due to the flux limited nature of the ASM sample and in- 
completeness of the optical identifications/distance mea- 
surements beyond ~ 10 kpc, the apparent luminosity func- 
tion which can be derived straightforwardly from the ASM 
flux measurements and the source distances (thin line his- 
tograms in Fig. [l2] and |l3|) needs to be corrected for the 
fraction of the Galaxy observable by ASM. This correction 
can be performed using the model of the spatial distribu- 
tion of X-ray binaries constructed in the previous section: 



dN 
~dL 



dN 

H 



M(< D(L)) 
Mua 



(11) 



Fig. 11. Fraction of the mass of the Galaxy visible to ASM 
with account for the selection criteria described in the text 
as a function of source luminosity. 



where ^ is the true luminosity function, (^x") obs - appar- 
ent luminosity function constructed using ASM flux mea- 
surements and the source distances, M(< D) - mass of the 
Galaxy inside distance D from the Sun computed using 
the volume density distributions for HMXB and LMXB 
sources from the Sect. ||, M to t - total mass of the Galaxy, 
D(L) is defined by: 



D(L) = min 



( L 



(12) 



where Fu m is the limiting (minimum) flux and D max - 
the maximum distance from the Sun of the sources used 
for constructing the luminosity function. As discussed in 
the previous sections we accepted the following selection 
criteria: Fu m = 0.2 cnts s _1 w 6.4- 10~ n erg s _1 cm~ 2 , i.e. 
equal to the completeness flux limit of the ASM catalogue, 
and D max = 10 kpc - a completeness limit of distance 



the mass fraction 



measurements estimated in Sect. 4.6. 

Obviously, for a given flux limit Fu 
of the Galaxy M ^ff[ L ^ is a decreasing function of the 
source luminosity as shown in Fig. |ll|. For the ASM sen- 
sitivity/completeness limit of w 6.4 ■ 10 -11 erg s _1 cm~ 2 
the entire volume inside D max = 10 kpc from the Sun is 
observable down to a luminosity of w 10 36 erg s" 1 (the 
flat part of the curves in Fig. [ll]) below which the mass 
fraction of the observable part of the Galaxy begins to de- 
crease. As the spatial distributions of HMXB and LMXB 
sources differ significantly, the volume correction and the 
luminosity function were calculated separately for HMXBs 
and LMXBs. The apparent and volume corrected (true) 
cumulative luminosity functions are presented in Fig. [L2[ 
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Fig. 12. The apparent (thin histogram) and volume corrected (thick histogram) cumulative luminosity function for 
LMXBs and HMXBs. The solid lines are the best fits to the data. 
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Fig. 13. The apparent (thin histogram) and volume corrected (thick histogram) differential luminosity function for 
LMXBs and HMXBs binned into bins with logarithmic with of 0.5. The solid lines are the best fits to the cumulative 
distributions. The fall-over of the apparent distributions below ~ 10 36 erg s _1 are due to the flux limited nature of 
the ASM sample (see Fig. [n]) 



Fig. [13| shows the corresponding differential distributions 
binned logarithmically over luminosity. 

The cumulative luminosity function of HMXBs (Fig. 
|l2| , right panel) does not seem to contradict to a power law 
distribution down to a luminosity of ~ 2 • 10 35 erg s _1 with 



some indication of flattening at lower luminosity. However, 
limited sensitivity of ASM and correspondingly large val- 
ues of the correction factor (Fig. [ll]) at low luminosities 
do not allow one to draw a definite conclusion regarding 
the shape of the luminosity function at these low lumi- 
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nosities (see comparison with ASCA source counts in Sec. 
0) . We therefore fitted the luminosity function of HMXBs 
in the L > 2 ■ 10 35 erg s -1 range with a power law distri- 



Model B: The same as the model A, except that 
the inner cut-off of the disk was set to r m = 4 kpc 
in accordance with the result of pchncn fc Binncy 



bution Using a Maximum-Likelihood method the best fit (1998)(dotted histogram in Fig. 14) 



parameters are: 



N(> L) =20- ( 



-0.64±0.15 



10 36 erg s 



(13) 



where L is the source luminosity in erg s 1 and N(> L) - 
total number of sources on the sky with luminosity greater 



Model C: The spheroid component is the same as in 
Model A. The disk radial distribution is without the 
inner cut-off, i.e. r m — and without modulation by 
the spiral structure. No bulge component is included 
for either LMXBs or HMXBs. The resulting density 



than It. 



distribution is similar to that derived by van Paradijs 



The shape uf buLh ciiiimlaLive and diffeieiilial lumiiius- 



ity function for LMXBs (Figs. [L2| |13|, left panels) indicates 
the presence of a high luminosity cut-off. We fitted the un- 
binned cumulative distribution with the functional form 



N(> L) = A - (L- a - L 



-a \ 
id ax J 



(14) 



corresponding to a power law differential luminosity func- 
tion with a sharp cut-off at L rnax . The value of the cutoff 
was set equal to to 2.7 • 10 38 erg s _1 which corresponds 
to the luminosity of the most luminous source within 10 
kpc, Sco X-l. The best fit values of other parameters are: 



N{> L) = 105- (( 



10 36 erg s' 



-) 



-0.26±0.08 _ 27Q-0.26\ 



(15) 



Note that the smaller number of sources and the 
steeper slope of luminosity function make the HMXB data 
insensitive to a high luminosity cut-off above ~ few x 10 36 
erg s -1 . 

5.1. Effect of the Galaxy model on the luminosity 
function 

From Eq. (|ll]) it is clear that the luminosity function de- 
pends on the spatial distribution of XRBs in the Galaxy. 
As discussed above, using the distance measurements 
available, we were able to determine some of the param- 
eters of their distribution. But the data are not sufficient 
to determine the entire distribution unambiguously. Thus 
we had to assume a spatial distribution of XRBs in the 
Galaxy. In order to investigate the effect of the adopted 
spatial distribution of X-ray sources on the derived lumi- 
nosity function we varied our model and computed the 
respective luminosity functions. 

For our analysis we used three different distributions 
for LMXBs and HMXBs. In the case of HMXBs, only the 
disk component was included in each of the three distri- 
butions. The modulation of the disk distribution by the 
spiral pattern, when present, was 100% for HMXB and 
20% for LMXB. The models are: 

— Model A: Our primary model constructed in Sec. 
and used to derive the luminosity function above 



(shown as a solid histogram in Fig. 14). 



fc White (1995[ ) for NS LMXBs (dashed histogram in 
Fig.|I4|). 



The resulting luminosity function for each of the three 
models are shown in Fig. [l4|. It is obvious that there is 
no strong dependence of the luminosity function on the 
mass distribution. The slopes vary in the range from 1.28 
- 1.30 for LMXBs and 1.64 - 1.72 for HMXBs. The total 
number of sources varies from 88 to 90 for LMXBs and 
from 21 to 26 for HMXBs. It is worth noting that the 
spiral pattern is no significant factor in the determination 
of the luminosity function of HMXBs although the spatial 
distribution shows clear signs of them. 

5.2. Total X-ray luminosity of Galactic X-ray binaries 

The total luminosity of all X-ray binaries in the Galaxy 
is calculated in the following way. Down to a luminosity 
of 10 36 erg s _1 we sum the measured luminosities of the 
individual sources to obtain a more precise number. For 
the lower luminosities that contribute only a small fraction 
to the total luminosity we use the analytical description 
of the luminosity function given by Eqs. (|l|) and @. 

The integrated luminosity of HMXBs and LMXBs in 
the 2-10 keV ASM band calculated in such way are ~ 



2-10 38 erg s _1 and: 



2.5 ■ 10 39 ergs s 1 , respectively. Note 



that these numbers refer to the luminosity averaged over 
the period from 1996-2000. The variability of individual 
sources or an outburst of a bright transient can change 
the luminosity by a factor of up to ~ 2 — 3. Due to the 
shallow slopes of the luminosity functions the integrated 
X-ray emission of the Milky Way is dominated by the 
~ 5 — 10 most luminous sources (see Table | and |). The 
maximum and minimum values for the luminosities were 
estimated by eye from the 1 day averaged light curves. 
The values in the tables therefore differ from the values in 
Table |. 

Normalised to the star formation rate which is about 



4 M Q yr" 1 in the Milky Way ( |McKee fc Williams 1997| ) 



galactic HMXBs emit about ~ 5 • 10 37 erg s _1 /(M Q yr" 1 ) 
The luminosity of LMXBs normalised to the stellar mass 
is about ~ 5 • 10 28 erg s^ 1 Mg 1 , assuming a stellar mass of 
the Galaxy of about 5 • 10 10 M Q . 

The contribution of Be X-ray binaries from the ASM 
sample to the integrated luminosity of HMXBs is ~ 5%. 
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Fig. 14. Dependence of the luminosity function on the adopted model of the spatial distribution of XRBs. The figures 
show the luminosity functions of LMXBs (left panel) and HMXBs (right panel) for three different Galaxy models. The 
solid, dotted and dashed lines in both panels correspond to the models A, B and C. 



Note that poor knowledge of the shape of the luminos- 
ity function at low luminosities, L & 10 35 erg s _1 should 
not influence the total luminosity considerably unless the 
luminosity function steepens significantly at these low lu- 
minosities (see Sec. [?]). 

The total number of X-ray binaries above 2 • 10 35 erg 
s _1 obtained from the luminosity functions is about ~ 190 
of which - 55 are HMXBs and ~ 135 - LMXBs. 



5.3. Luminosity function and M distribution of X-ray 
binaries 

The X-ray luminosity function is obviously related to the 
distribution of X-ray binary systems over the mass loss 
rate of the secondary, M . The simplest assumption would 
be that both distributions have the same slope in the 
range corresponding to luminosities of ~ (0.01 — l)LEdd- 
At larger luminosities, L ^ LEdd, the luminosity function 
has a break or cut-off, well in accordance with theoretical 
expectation, that the luminosity due to accretion cannot 
exceed the Eddington luminosity of the primary star by 
a large factor (see discussion in Sect. ||). The donor star 
in a binary system, on the other hand, "does not know" 
about the Eddington critical luminosity, therefore the dis- 
tribution of binary systems over the mass loss rate of the 
secondary, M, is not expected to break near the Eddington 
value for the compact object. Thus the distribution of bi- 
nary systems over M is expected to continue with the 
same slope well beyond the Eddington value. 

Extremely super-Eddington values of the mass accre- 
tion rate M can result in quenching of the X-ray source 



and/or its obscuration by the matter expelled from the 
system by radiation pressure (Shakura & Sunyaev 1973). 
This would lead to the appearance of a peculiar object, 
dim in X-rays and extremely bright in the optical and 
UV band - similar to SS 433 or the recent fast transient 
V4641 Sgr at the peak of its optical outburst. Such ob- 
jects would emit only a negligible fraction in the X-ray 
band and would contribute to the lower luminosity end of 
the XRB luminosity function. 

For moderately super-Eddington values of M 10 — 
lOOMEdd, however, one might expect the appearance of 
a near- or slightly super-Eddington source, therefore all 
such systems are expected to cluster near LEdd- For a 
given slope of the luminosity function the number of such 
sources can be easily estimated. For the observed param- 
eters of the LMXB luminosity function (slope = 1.3, 42 
sources with 36.5 < \og(Lx) < 38) and assuming that 
the M distribution continues with the same slope = 1.3, 
the total number of sources with M corresponding to the 
range of luminosities of 10 38 — 10 39 and 10 39 — 10 40 erg 
s _1 is w 10 and « 6 correspondingly (w 7 sources are 
expected to have M corresponding to L > 10 40 erg s _1 ). 
These estimates are in disagreement with the actually ob- 
served number of sources with L <; 10 38 erg s , which is 
equal to 8. In order to reconcile the expected number of 
sources near LEdd with the observations, a slope of the M 
distribution of <; 1.35 — 1.40 is required which is somewhat 
steeper than the observed value of ~ 1.3. We note that the 
slope of ~ 1.35 is within ~ lcr of the the observed value. 

Finally, there are several effects that can suppress the 
number of the low luminosity sources, i.e. make the lumi- 
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Table 5. List of the most luminous LMXB sources con- 
tributing w 90% to the integrated luminosity of LMXBs 
in the 2-10 keV band, averaged over 1996-2000. The 12 
most luminous sources contribute « 80% of the integrated 
luminosity. 



Table 6. List of the most luminous HMXB sources that 
contribute ss 40% to the integrated luminosity of HMXBs 
in the 2-10 keV band, averaged over 1996-2000. 



Source 



Cir X-l 

GRS 1915+105 
Sco X-l 
Cyg X-2 
GX 349+2 
GX 17+2 

GX 5-1 
GX 340+0 
GX 9+1 
NGC 6624 
Ser X-l 
GX 13+1 



X 1735-444 



XTE J1550-564 
KS 1731-260 
X 1705-440 



X 1624-490 



L x [10 JB erg s" 1 ] 



avg. 



4.4 
3.7 
2.7 
1.8 
1.6 
1.5 

1.4 
1.3 
0.75 
0.47 
0.43 
0.41 



0.3 

1 
2 

0.9 
1.1 
1.1 

1 

0.9 
0.5 
0.15 
0.26 
0.25 



10 
11 
4.5 
3.4 
2.7 
2.4 

1.8 
1.8 
1.0 
0.8 
0.6 
0.6 



0.35 



0.2 



0.6 



0.35 
0.28 
0.25 



0.005 
0.06 
<0.04 



2.1 
0.6 
0.6 



0.24 <0.13 



0.4 



dist. 

[kpc] 



10.9 
12.5 
2.8 
11.3 
9.2 
9.5 

7.2 
11.0 
7.2 
8.0 
8.4 
7.0 



9.2 



5.3 
8.5 
7.4 



13.5 



Ref. 



1 

2 
3 

1,4-7 

1,7,8 

1,7,9, 

10 

1,7 

1,7 

9 

10,11 

12 

13 



14 

15-17 
18 







min. and max. luminosity were estimated by eye from the 
1 day averaged light curves. 
References for the distances: 
(2) 



Source 


L x 


fl0 38 ere 


° J 


dist. 


Ref. 




avg. 


min.M 


max (") 


[knrl 






0.5 


0.08 


1.4 


9.0 


1 


Cpn X-^ 


0.15 


<0.03 


0.7 


9.0 


2-5 


Cyg X-l 


0.05 


0.02 


0.17 


2.1 


6 


X 1657-415 


0.043 


<0.02 


0.22 


11.0 


7 


V4641 Sgr 


0.028 


<0.02 


7.3 


9.9 


8 


GX 301-2 


0.02 


< 0.005 


0.4 


5.3 


9 


XTE J 1855-024 


0.015 


<0.01 


0.11 


10.0 


10 


X1538-522 


0.014 


< 0.008 


0.08 


6.4 


11 


GS1843+009 


0.01 


< 0.007 


0.11 


10.0 


12 


X1908+075 


0.008 


< 0.006 


0.05 


6.4 


13, 
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min. and max. luminosity were estimated by eye from the 

1 day averaged light curves. 

References for the distan ces: (1) - [Predehl ct al. (200C ) , (2 ) - 
[Krzeminski (1974| ) , (3) -|Hutchings et al. (1979D, (4) 



et al. 



(1997|), (5) 



Bahcall (197: 



(6) - [Massey et al. (1995 ) 
(7) - phakrabarty~eT al. (1993Q , (8) - Prosz et al. (2000 ) 
(9) -|Kaper et al. (1995[ ), (10 ) - porbet et al. (1999Q , (11 



Motch 



|Reynolds et al. (1992|)~(12) - [Israel et al 
et al. (200q ), (14) - |van Paradijs fc White 



+99 



2001|), (13) 



Wen 



( 4 a 



( 1) -|van Paradijs fc White ( 1995 ) 



Smale (1998) 



(9) ■ 
OIL 



Webbink (1985|), (12) - [Ebisuzaki et al. (1984|) 



(1996 h. 
(1993^ 

(13) - 
et al. (] 

(18) -paberl fc Titarchuk (1995|) 



(7) -|Pennirix(1989|), (8) - [Wachter fc Marg 
Christian fc Swank (199%, (10) - pjorgovsk 



Vlirabcl fc Rodriguez (1994| ), (3) -[Bradshawet al. (1999j ), 
Orosz fc Kuulkers ( 19991)7(5) - |Co"wley et al. (197S|), (6) 



Bandyopadh yay et al. (1999|), (14) - ?, (15) - |Barret | 
al. (1997|), (17) - punyaev (lOOCf ) 



998|). (16) -pmith et 



nosity function flatter than the M distribution. The most 
obvious and important are discussed below. 

— In the case of HMXBs the magnetosphere of the 
strongly magnetised, rapidly rotating neutron star can 
prevent the accretion at low M via the propeller effect 
( Illarionov fc Sunyaev 1975 ). 

— Be-systems are characterised by regular outbursts cor- 
responding to the passage of the neutron star through 
the equatorial stellar wind. Therefore for such sources 
the true value of the M in the binary system is mea- 
sured by the peak luminosity during the outbursts 
whereas the long term averaged luminosity, used to 
construct the luminosity function, can give a signifi- 
cantly underestimated value. 

— A common property of LMXBs, containing both neu- 
tron stars and black hole, is the presence of relativis- 
tic jets which might carry away a sizable fraction of 
the energy of accretion ( Mirabel fc Rodriguez 19*9^ ). 
The presence of jets correlates with the X-ray spectral 
state: the jets are absent (and hence the true accretion 
efficiency is higher) in the soft spectral state corre- 



sponding to higher valu es of M. The jets exist only in 
the hard spectral state ( Fender 2001 ), thus decreasing 
the accretion efficiency at lower M. 

— In the case of black hole binaries an ADAF can form at 
low accretion rate in which case the accretion efficiency 
is proportional to M and the X-ray luminosity scales 
as L oc M 2 ( |Narayan fc Yi 1995| ). 

— At sufficiently low accretion rates a source becomes a 
transient with a recurrence time varying from ~ 1 to 

50 years (White et al. 1984). This would decrease 
the number of low and intermediate luminosity sources 
in the luminosity function constructed on the several 
years baseline. 

The number of luminous X-ray binaries in the Milky 
Way is insufficient to study the shape of the luminosity 
function near LEdd in detail. On the other hand within 
next several years CHANDRA X-ray observatory will 
study compact sources in a large number of nearby, d Si 50 
Mpc galaxies and the total number of the X-ray binaries 
detected in other galaxies can easily reach several hundred 
or thousand. In this context it might be interesting to con- 
struct a combined luminosity function of X-ray binaries in 
our and other galaxies to study its exact shape at the high 
luminosity end and search for a possible excess of sources 
near L E dd- 

6. Comparison with nearby galaxies 

The total luminosity of X-ray binaries in the Milky Way, ~ 
2 — 3- 10 39 erg s _1 in the 2-10 keV band, agrees sufficiently 
well with observations of M 31, for which GINGA has 
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Fig. 16. Cumulative luminosity functions of galaxies observed with CHANDRA. The left panel shows actively star 
forming spiral galaxies that include NGC 4038/39 and M 82 which are supposed to be dominated by HMXBs. For 
comparison the luminosity functions of Galactic X-ray binaries and HMXBs alone are shown. The right panel shows 
elliptical galaxies including the SO galaxy NGC 1553. For comparison the luminosity function of Galactic LMXBs is 
shown. 
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Fig. 15. 

LMXBs 



Cumulative luminosity function of Galactic 
and also the best fit values for the XMM-Newton 



observ ation of M 31 by SMrgy et al. (2001 ) 



Recently XMM-Newton observed the inner 30' region 
of M 31 (Shirey et al. 2001). In total 116 sources were de- 
tected above the limiting luminosity of 6 • 10 35 erg s _1 in 
the 0.3-12 keV energy range, assuming a distance of 760 
kpc. Shirey et al. (2001) distinguish between two luminos- 
ity ranges, 36.2 < log(Lx) < 37.4, for which the best fit 
slope is -0.47 ± 0.03, and 37.4 < log{L x ) < 38.1 where 
the best fit slope is — 1.79±0.26. At the distance of 760 kpc 
30' correspond to ~ 6.6 kpc therefore these data should 
be compared with the luminosity function of Galactic 
LMXBs, assuming that similarly to the Milky Way the 
inner part of M31 is populated mainly with LMXBs. The 
two luminosity functions are plotted in Fig. Although 
the general shapes of the luminosity functions of LMXBs 
in the Milky Way and in M 31 are similar, it is obvious 
that one can not be obtained from the other by a shift 
along the vertical axis as one would expect if the luminos- 
ity function was simply proportional to the mass of the 
host galaxy. 

CHANDRA observations have produced luminosity 
functions of compact sources in a number of nearby galax- 
ies, in cluding ellipticals: NGC 4 697 ( |Sarazin et al. 2000|) 
M 84 (|Finoguenov fc Jones 200l|) and NGC 1553 ([Blantor 



et al. 2001), spirals: M 81 (Tennant et al. 2001), Circinus 



( |Smith fc Wilson 200l[ ), M31 ( parcia et al. 2000[) and 
starburst galaxies: NGC 4038/39 (Antennae) (Fabbianc 



et al. 2001) and M 82 (Griffiths et al. 2000). The lumi 



found a luminosity of 5 ■ I0 iy ergs s 1 between 2-20 keV nosity functions of the compact sources in these galaxies 

are compared to that of the Milky Way in Fig. [l^. The 



(Makishima et al. 1989) 
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left panel in Fig. shows spirals and starbursts which 
are expected to have a higher fraction of HMXBs due to 
higher star formation rates. These are compared with the 



In order to study the low flux regime below the ASM 
completeness limit of as 6.4 • 10~ n erg s _1 cm~ 2 , we use 



ASCA data from the Galactic Ridge Survey (Sugizak 



luminosity functions of HMXBs and all X-ray binaries in 
the Milky Way. The right panel in Fig. |16| shows ellipti- 
cal galaxies along with the luminosity function of Galactic 
LMXBs. 

As the example of our Milky Way shows, X-ray binaries 
in globular clusters play an important role in determina- 
tion and understanding the properties of the population. It 
is also well known that globular cluster systems are quite 
different for early- and late-type galaxies, in terms of num- 
ber per galaxy luminosity ( Harris fc Racine 197S| ) as well 



as depend on the environment of the host galaxy (Bridges 
& Han ;s 1990| ). Taken together this shows the need for a 
closer study of X-ray binaries in globular clusters - ide- 



et al. 2001) covering ss 40 square degrees with the lim- 
iting sensitivity of ~ 3 • 10~ 13 erg s _1 cm~ 2 . Since most 
of the sources in the ASCA survey are unidentified we fol- 
lowed the criterion suggested by Sugizaki et al. (2001[ ) in 
order to discriminate X-ray binary candidates from other 
sources: that X-ray binary candidates have either a spec- 
tral photon index T < 1, or a spectral photon index T < 3 
and a column density Njj < 0.8 • 10 22 cm~ 2 . Excluding 
otherwise identified sources with these spectral properties 
there remain 28 sources. We fit the Log(N)-Log(S) of the 
selected sources with the procedure similar to that used 
for ASM sources, modified to account for the flux depen- 



dent sky cov erage of the ASCA survey (Fig. 7 in [Sugizak 



ally th fy should be treated separately, when studying the 
luminosity function of LMXB sources. Unfortunately only 
for few galaxies there are observations which allow the 
separation of globular cluster X-ray sources, e.g. M 31 ([D] 
Stefant j) et al. 2002[ ) and NGC 1399 ( |Angclini et al. 2001| ). 
We therefore decided to ignore in the present study the 
possible effects of the globular cluster sources on the over- 
all luminosity function. 

Comparing the HMXB luminosity function in our and 
nearby star forming galaxies we could check the propor- 
tionality of the HMXB luminosity to star forming rate. 
There might be several additional factors involved includ- 
ing chemical abundance of the particular galaxy. For ex- 
ample, the HMXB sources in LMC and SMC appear to 
be significantly more luminous than the HMXB sources in 
our Galaxy, even though the star formation rates are com- 
parable. Especially interesting is the case of the Antennae 
galaxies where the difference from the Galactic HMXB lu- 
minosity function is extremely impressive. It seems that it 
can not be explained simply by the difference in the star 
format ion rate , which is about 20 times higher ( Neff fc 
Ulvesti id 2000 ) whereas the number of X-ray sources is a 
factor of more than 50 higher. This example shows that 
the knowledge of the HMXB luminosity function seems to 
be insufficient to measure the star formation rate in galax- 
ies and to estimate the distances to them with acceptable 
precision. 

CHANDRA observations are also opening an impor- 
tant possibility to check the proportionality of LMXB lu- 
minosity functions to the mass of the parent galaxies. 



7. Low luminosity sources 

7.1. Extension of Log(N)-Log(S) towards lower fluxes 

Since the sensitivity of ASM is limited to relatively high 
flux sources it is interesting to investigate the behaviour 
of the Log(N)-Log(S) at lower fluxes. Note that, given 
the slope observed by ASM (1.2 and 1.61 for LMXBs and 
HMXBs), the Log(N)-Log(S) distribution should flatten 
at low fluxes since the total number of sources in the 
Galaxy is finite. 



et al. (20011) ). The resulting Log(N)-Log(S) is: 



N(> S) = 9.4-10" 



-0.42±0.08 



(16) 



where S is flux in units of 10 -12 erg s _1 cm~ 2 . To compare 
ASCA data with an extrapolation of the ASM number- 
flux relation one needs to account for the difference in 
their sky coverage (|Z| <, 40° and \b\ ^ 0.3° for ASCA 
survey and entire sky for the ASM data). An approxi- 
mate value of the correction factor can be estimated as 
the fraction of the mass of the Milky Way covered by the 
ASCA Galactic Ridge Survey with account for its sensitiv- 
ity and the particular pattern of its sky coverage (Fig. 1 in 
|Sugizaki et al. (200l| )). The mass fraction was calculated 
using the Galaxy model described in Sect. || and equals 
to ~ 1:21. Converting the cumulative Log(N)-Log(S) to 
differential Log(N)-Log(S) for ASCA X-ray binary can- 
didates and all ASM X-ray binaries and multiplying the 
resulting ASCA Log(N)-Log(S) by 21 we obtain the result 
shown in Fig. |l7|. 

It is obvious that the agreement between ASM and 
ASCA data is sufficiently good. The slopes are different at 
the ~ 2<T level. On the other hand since the sources are all 
unidentified and their distances unknown it is not possible 
to distinguish between high and low mass X-ray binaries 
which have different slopes of their Log(N)-Log(S) dis- 
tributions in the ASM sample. Indeed, due to the small 
range in Galactic latitude b 11 covered by the ASCA survey 
and due to the fact that HMXBs have a 3 times smaller 
vertical scale height (cf. Sec. ||), the ratio of HMXBs to 
LMXBs should be different for the ASCA and ASM sam- 
ples. The fraction of HMXBs, having steeper Log(N)- 
Log(S), should be larger in the ASCA sample and thus 
the resulting Log(N)-Log(S) should be somewhat steeper. 
We conclude that the data of the ASCA Galactic Ridge 
Survey indicate that there are no significant deviations in 
the Log(N)-Log(S) from the extrapolations of the ASM 
data down to the sensitivity limit of the ASCA survey of 
erg s _1 cm" 
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7.2. Low luminosity end of X-ray binary luminosity 
function 

Knowledge of the Log(N)-Log(S) observed by ASCA and 
the spatial distribution of sources in the Galaxy gives a 
possibility to constrain the low luminosity end of the lumi- 
nosity function. If the luminosity function observed with 
ASM continues to lower luminosities then it should be pos- 
sible to reproduce the Log(N)-Log(S) observed by ASCA 
according to the formula 



N(> S) 



dN M(< r) 



ASCA 



dL 



M, 



total 



with 



dL, (17) 



(18) 



where iV(> S) is the number of sources with a flux higher 
than S observed by ASCA, 4^ is the differential luminos- 



M(<i-)asca 



is the fraction of mass within 



ity function, and , ( 
a radius r from the Earth within the field of view of the 
ASCA survey, L max is the high luminosity cut-off of the 
luminosity function (Eqs. ( |l3| ) and (|l5|)). The L min is the 
low luminosity cut-off of the luminosity function below 
which it is assumed to be equal to zero. This quantity 
characterises roughly the luminosity level at which the lu- 
minosity function deviates significantly from the extrapo- 
lation of the ASM power law. 
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Fig. 17. Comparison of the differential Log(N)-Log(S) re- 
lation for Galactic X-ray binaries obtained by ASM (solid 
line with break) and by ASCA Galactic Ridge Survey 
(dashed line). The ASCA number-flux relation was mul- 
tiplied by an approximate correction factor accounting for 
the difference in the sky coverage of the ASM and ASCA 
surveys (see text for details). 
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Fig. 18. Comparison of the number-flux relation observed 
in the ASCA Galactic Ridge Survey (points) and the pre- 
dicted number-flux relation based on the extrapolation of 
the ASM luminosity function to low luminosities (lines). 
The vertical axis shows the number of sources in the en- 
tire field of the ASCA survey. The ASCA number-flux 
relation was corrected for the flux dependent sky coverage 
(Fig. 7 in Sugizaki et al. (2001 )). The predicted number- 
flux relations were computed according to Eq. (^) using 
the extrapolation of the ASM luminosity functions and the 
volume density distributions of X-ray binaries described in 
Sect. [|. The thick solid lines show the combined Log(N)- 
Log(S) of LMXBs and HMXBs for different values of the 
low luminosity cut-off. The thin dashed lines show the con- 
tributions of LMXBs and HMXBs separately for the case 
without cut-off. 



The predicted Log(N)-Log(S) calculated from Eq. (|17|) 
is compared with the Log(N)-Log(S) of X-ray binary can- 
didates from the ASCA survey in Fig. [l^. In plotting the 
ASCA data (solid circles) we added five bright sources lo- 
cated in the ASCA field of view that were excluded from 



the final catalogue in Sugizaki et al. (2001 ) and corrected 
for the flu x dependent sky cov erage of the ASCA survey 
(Fig. 7 in |Sugizaki et al. (200l| )). The predicted Log(N)- 
Log(S) was calculated according to Eq. ( p~7| ) separately for 
HMXB and LMXB using the extrapolation of the respec- 
tive ASM luminosity functions. The mass integral M(< r) 
in Eq. (^) was calculated taking approximately into ac- 
count the actual pattern of ASCA pointings and using the 
volume density distributions constructed in Sect. ^. The 
predicted combined Log(N)-Log(S) of HMXB and LMXB 
sources is shown in Fig. [l8] by the thick solid lines for dif- 
ferent values of the low luminosity cut-off L m i n - The thin 
solid and dashed lines show the contributions of HMXBs 
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and LMXBs respectively for the case without low lumi- 
nosity cut-off. 

It is clear from Fig. [IJ that the predicted number- 
flux relation of X-ray binaries agrees with the ASCA data 
very well. Given the volume density distributions of X- 
ray binaries in the Galaxy, the low flux end of the ASCA 
Log(N)-Log(S) is sensitive to sources with luminosities of 
~ 10 34 erg s _1 . The good agreement with the predicted 
Log(N)-Log(S) distribution implies that the data do not 
require a low luminosity cut-off of the luminosity function 
down to ~ I0 34 erg s _1 . 

7.3. Young objects in star forming regions 

Recent observations with the CHANDRA X-ray observa- 
tory of the Orion Nebula cluster allow one to estimate the 
contribution to the X-ray emission from young objects in 
the star forming regions. |Schulz et al. (200l| ) observed the 
Orion Trapezium region and found 111 sources above the 
sensitivity threshold of 6.6 • 10 28 erg s^ 1 , assuming a dis- 
tance of 440 pc. The total luminosity of their sample is 
about 5.6 • 10 32 erg s _1 . This luminosity is dominated by 
the brightest source in the Orion Nebula cluster, 8 1 Ori 
C, which provides about 1.8 • I0 32 erg s _1 . Extrapolating 
this result to the whole Orion Nebula Cluster in which 
CHANDRA observed about 1000 sources we obtain a to- 
tal luminosity of the star cluster of about 4 • 10 33 erg s _1 , 
counting the luminosity of 8 1 Ori C only once and mul- 
tiplying the rest by 10, assuming the luminosity function 
of the Trapezium region is representative for the whole 
Orion Nebula cluster. To estimate the X-ray luminosity 
of all star forming regions in the Galaxy one can proceed 
in two ways. Taking the mass of the molecular gas in the 
Orion cluster to be ~ 1O 5 M ( [Maddalena et al. 1986] ), 
and the total mass of the molecular gas in the Galaxy 



to be ~ 10 9 M o (Williams fc McKee 1997), the total lu- 



y-ts- 



4 ■ 10 



org c 



On the other hand one 



mmosi 

can us h the star formation rate in the Orion Nebula clue 
ter and the Galaxy as the determining factor. Taking the 

O yr 



SFR in Orion to be > 10~ 4 M G) yr" 1 ([Hillenbrand 1997|) 



and the SFR in the Galaxy to be 4 M yr _T ([McKee fe 



Williams 1997), the total luminosity of young objects in 



the star forming regions in the Galaxy is 1.6 • 10 erg 
s _1 . Taking into account that the latter value is an up- 
per limit, both numbers agree sufficiently well. Therefore 
star forming regions contribute less than ~ few per cent 
to the integrated X-ray emission of the Galaxy but ~ 20% 
or more to the luminosity of HMXBs in the energy range 
from 2-10 keV. On the other hand the spectrum of young 
stellar objects is much softer than the spectrum of X-ray 
binaries. 



8. High luminosity sources 

In recent months the CHANDRA X-ray observatory was 
able to resolve single X-ray sources in other galaxies that 
appear to radiate at or above the Eddington limit for a 1.4 
M neutron star, i.e. ~ 2 • 10 38 erg s _1 . Similar behaviour 



is also observed in Galactic X-ray binaries by ASM. The 
slightly different spectral band used in these CHANDRA 
observations, usually 0.3-10 keV compared to 2-10 keV for 
ASM, does not lead to significant differences in luminosity. 

Table [l] lists the sources which were observed either by 
ASM or some other instrument to emit at or above the 
Eddington limit for a 1.4 M Q neutron star. The spatial 
distribution of these sources is shown in Fig. [l9| and can 
be compared to the distributions of the brightest sources 
observed by CHANDRA in other galaxies. 

There are several reasons why sources can emit super- 
Eddington luminosity: 

— For accreting black holes in high state radiation is com- 
ing from the quasi-flat accretion disk where electron 
scattering gives the main contribution to the opacity. 
Under these conditions the radiation is emitted accord- 
ing to 



f( P ) = (1 + 2.08 M )/x 



(19) 



where /i — cos(i) where i is the inclination angle. It is 
easy to show that the radiation flux perpendicular to 
the plane of the disk exceeds the average value by 3 
times (see Shakura fc Sunyaev (1973| ) for discussion). 
Some of the normal stars entering the X-ray binary 
phase are strongly evolved and have an unusual chem- 
ical abundance, e.g. if a He-enriched star supplies mat- 
ter the Eddington luminosity is twice higher than for 
hydrogen plasma due to the change in cross-section per 
nucleus. 

Just these two factors permit to surpass the classical 
Eddington limit by a factor of ~6. 
The star supplying material to the neutron star or 
black hole "does not know" about the existence of the 
Eddington luminosity limit due to accretion. Therefore 
some part of the matter will outf low forming a super- 
critica l disk. In the approach of Shakura fe Sunyaev 
(1973 ) it is possible to gain a factor of ln(m) ss 3 — 5 for 

M 



>> 1 with m = 



Paczynsky fc Wiita (1980| ) 



and Abramowicz et al. (1 



constructed the solution 
of slim disks which also permits luminosities higher 
than the Eddington luminosity. 

Many X-ray binaries show from time to time the accel- 
eration of powerful jets ( Mirabel fc Rodriguez 1999 ). 
These relativistic jets might produce strongly beamed 
X-ray emission with flux strongly exceeding the aver- 
age and Eddington critical value for isotropic sources. 
See also the discussion by Kocrding et al. (2001 ), 



Fabrika fc Mescheryakov (2000| ) and [King et al. (200l| ) . 
In the case of accretion on to a neutron star with 
strong magnetic field the accretion columns form near 
the surface of the neutron star in the polar regions. 
Such columns can have a super-Eddington luminosity, 
because photons are emitted perpendicular to the axis 
of the accretion column and the light pressure force is 



balanced by magnetic field (Basko & Sunyaev 1976). 
In Z-sources (luminous accreting neutron stars with 
low magnetic field) the boundary layer width expands 
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Fig. 19. The spatial distribution of Galactic X-ray binaries that have shown episodes of Eddington or super-Eddington 
luminosity for a 1.4 M Q neutron star. The coordinate system is the same as in Fig. [?]. Filled circles indicate HMXBs, 
open squares indicate LMXBs. Note that fact that the majority of the sources arc located at y > reflects the flux 
limited nature of the ASM sample. 



rapidly with increasing accretion rate reaching several 
star radii (Popham & Sunyaev 2001). This quasi-flat 
continuation of accretion disk might also have super- 
Eddington luminosity of the type of the slim disk. 

9. Summary 

We studied the population of X-ray binaries in the Milky 
Way. 



well described by a simple power law with a slope of 



earlier results ( 


van Paradijs & White 1995; 


White & 


van Paradijs 1996]; 


Koyama et al. 199C; 


Nagase 1989) 



we found significant differences in the spatial (3-D) dis- 
tribution of high and low mass X-ray binaries. HMXBs 
are more concentrated towards the Galactic Plane with 
a vertical scale height of 150 pc, tend to avoid the 
Galactic Bulge and central ~ 3 — 4 kpc of the Galaxy 
and show clear signatures of the spiral structure. The 
distribution of LMXB sources, on the contrary, peaks 
strongly at the Galactic Bulge and shows a pronounced 
minimum at ~ 3 — 4 kpc. Some signatures of the 
Galactic spiral structure are also present. The vertical 
distribution of LMXB sources is significantly broader, 
with a scale height of 410 pc. 

We constructed the long-term averaged Log(N)- 
Log(S) distribution of high and low mass X-ray bi- 
naries in the 2-10 keV energy range using the data of 
the ASM instrument aboard RXTE from 1996-2000 
to the limiting sensitivity of w 6.4 • 10~ n erg s _1 
cm" 2 . The Log(N)-Log(S) distribution of HMXBs is 



the differential distribution of 1.61 + o!i2 down to a flux 
10~ n erg s _1 enr 2 



limit of s» 6.4 • 10 erg s cm . The differential 
Log(N)-Log(S) distribution of LMXBs has a slope of 
— 1.2 ± 0.06 and requires a high-flux cutoff at ~ 110 
ASM cnts s _1 , w 3.5 • 10 -8 erg s _1 cm -2 . A compar- 
ison with data of the ASCA Galactic Ridge Survey 
(Bugizaki et al. 2001) which covered ~ 40 square de- 
grees with ~ 100 times better sensitivity did not reveal 
any evidence of significant departures of the Log(N)- 
Log(S) from an extrapolation of the ASM data down 
to w 5 • 10~ 13 erg s" 1 cm~ 2 . 

Using the source distances available and assuming a 
model for the volume density distribution we con- 
structed luminosity functions for HMXBs and LMXBs 
in the 2-10 keV energy range. The sensitivity limit of 
the ASM catalogue allows one to study the XRB lumi- 
nosity functions down to a luminosity of ~ 2 • 10 35 erg 
s _1 . The differential luminosity functions can be de- 
scribed by a power law with slopes of 1.64 and 1.27 for 
HMXBs and LMXBs respectively. For LMXB sources 
a cut-off at ~ 2.7- 10 38 erg s" 1 is required. The HMXB 
data are insufficient to detect a high luminosity cut-off 
above ~ few x 10 36 erg s" 1 . A comparison with the 
data of ASCA Galactic Ridge Survey did not find evi- 
dence for significant departures from these power laws 
down to luminosities of ~ 10 34 erg s _1 . 
The complete catalogue of our sample of X- 

at http: / / www.mpa- 



ray binaries is available 
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garching.mpg.de/~grimm/. Properties of the brightest 
sources are summarised in Tables ||, ||, |l|. 

— The integrated luminosity of X-ray binaries in the 
Milky Way in the 2-10 keV band averaged over 1996- 
2000 is - 2 - 3 ■ 10 39 erg a' 1 to which LMXB sources 
contribute ~ 90%. Normalised to the stellar mass and 
the star formation rate the integrated luminosity of 
LMXBs (~ 2.5 ■ 10 39 erg s" 1 ) and HMXBs (~ 2 • 10 38 
erg s _1 ) correspond to ~ 5 • 10 28 erg s _1 Mg 1 and 
~ 5 • 10 37 erg s _1 /(M Q yr" 1 ), respectively. The total 
number of the X-ray binaries brighter than 2 • 10 35 erg 
s _1 is ~ 190 of which ~ 55 are high mass and ~ 135 
are low mass binaries. Extrapolating the luminosity 
functions towards low luminosities we estimate the to- 
tal number of the X-ray binaries brighter than 10 34 
erg s" 1 as ~ 705 (~ 325 LMXB and - 380 HMXB 
sources). These estimates might be subject to the un- 
certainty of a factor of ~ 2 due to insufficient knowl- 
edge of the spatial distribution of X-ray binaries in the 
Galaxy. 

— Due to the shallow slope of the luminosity function the 
integrated X-ray emission of the Milky Way is domi- 
nated by ~ 5 — 10 brightest sources. Variability of in- 
dividual sources or an outburst of a bright transient 
source can increase the integrated luminosity of the 
Milky Way by as much as a factor of ~ 2. 

— We found that at least 16 sources in the Galaxy showed 
episodes of super-Eddington luminosity for a 1.4 M Q 
neutron star. We plotted the distribution of these 
sources across the Galaxy in various projections, which 
can be used to compare with the recent CHANDRA 
and XMM-Newton images of the nearby galaxies. 
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